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Abstract 

Subspace codes were introduced by Kotter and Kschischang for error control in random linear 
network coding. In this paper, a layered type of subspace codes is considered, which can be viewed as 
a superposition of multiple component subspace codes. Exploiting the layered structure, we develop two 
decoding algorithms for these codes. The first algorithm operates by separately decoding each component 
code. The second algorithm is similar to the successive interference cancellation (SIC) algorithm for 
conventional superposition coding, and further permits an iterative version. We show that both algorithms 
decode not only deterministically up to but also probabilistically beyond the error-correction capability 
of the overall code. Finally we present possible applications of layered subspace codes in several network 
coding scenarios. 

Index Terms 

Network coding, error correction, subspace codes, superposition coding 

I. Introduction 

In the paradigm of network coding [H, information transmission is highly susceptible to packet errors. 
Due to the mixing nature, even a single corrupt packet may cause widespread error propagation, rendering 
the entire transmission useless. Thus error control is essential for providing reliable transmission in 
network coding. 

The notion of network error correction coding was first introduced by Cai and Yeung in t3]. Their 
approach is based on a coherent transmission model, in which both the transmitter and receiver know 
the network topology. 
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In the context of random linear network coding [2 ], Kotter and Kschischang proposed the subspace 
coding method as the error control solution |4|. A noncoherent transmission model was assumed where 
neither the transmitter nor receiver have knowledge of the network topology and the particular network 
codes used. Subspace codes encapsulate network codes to provide an end-to-end error protection. 

Recently, a coding scheme consisting of a number of subspace codes was proposed by Siavoshani et al. 
for multi-source multicast network coding Q. In (H, Dikaliotis et al. extended this work by constructing 
capacity-approaching subspace coding schemes for multi-source network coding transmission. 

In this paper, we investigate the superposition property of the codes in [7] and propose two decoding 
algorithms. Due to their layered structure, we refer to the codes as layered subspace codes. Our main 
contributions can be summarized as follows. 

• We provide more insights by showing that a layered subspace code forms a superposition coding 
scheme ifTUl . 

• We develop two efficient decoding algorithms. The first algorithm operates by separately decoding 
each component code. The second algorithm is similar to the successive interference cancellation 
(SIC) algorithm for conventional superposition coding, and further permits an iterative version. We 
show that both algorithms are guaranteed to decode up to the error-correction capability of the 
overall code. Besides, they can occasionally decode beyond the capability. 

• We point out that layered subspace codes can find more applications than presented in Q. For 
example, the codes can be used as an adaptive transmission scheme or an unequal error protection 
scheme for single-source multicast network coding. 

The rest of the paper is organized as follows. Section II gives a brief review of subspace codes. In 
Section III, we investigate the properties of layered subspace codes and develop two decoding algorithms 
for these codes. Section IV discusses some possible applications of layered subspace codes in network 
coding. Finally, we conclude the paper in Section V. 

II. Preliminaries 

In this section, we briefly recall the subspace coding method [4] for random linear network coding 
(RLNC) (21. In RLNC, a source injects some packets into the network, each being regarded as a row 
vector over a given finite field. These packets propagate though the network, passing though a number 
of intermediate nodes between source and receiver. Each intermediate node creates a random linear 
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combination of packets it received, and transmits this combination. Finally, a receiver collects a set of 
such randomly generated packets and tries to recover the packets injected into the network. 

Let F q be a finite field with q elements, where q is a prime power. Let W be a fixed finite-dimensional 
vector space over F g and V(W) the set of all subspaces of W. Denote by dim(V) the dimension of an 
element V G V(W). Two operations on V(W) can be defined ifTTTl . The intersection of V, U G V(W) is 
defined as 

V nu = {w.w gV,w eU}, (l) 

which is the subspace of largest dimension contained in both V and U. The sum of V and U is defined 

as 

V + U = {v + u : v G V,u G U}, (2) 

which is the subspace of smallest dimension containing both V and U. If V and U intersect trivially 
(i.e., V n U = {0}), V + U is called the direct sum, denoted by V U. 

For RLNC, the transmission is modeled as an operator channel, where both the input and output are 
a subspace of W (H. Let V be the input and U the output, the operator channel relates them by 

U = (V n U) © E, (3) 

where E is called the error space. In transforming from V to U, it is said that the operator channel 
commits p = dim(V) — dim(y n U) erasures (also called deletions of dimension) and t = dim(i?) = 
dim(f7) — dim(VnC7) errors (also called insertions of dimension). In practice, the source sends a basis for 
the information-carrying vector space V and the receiver collects a set of vectors that span the possibly 
corrupt vector space U. 

To measure the degree of dissimilarity between V and U, the subspace distance has been introduced 
[4] 

d s (V, U) = dim(V + U)- dim(V n U) 

= dim(F) + dim(C7) - 2dim(F n U) 

= p + t. (4) 
With the definition, V(W) forms a metric space. 
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A subspace code C is defined to be a nonempty subset of V(W) (4). Each codeword of C is a subspace 
of W. The minimum (subspace) distance of C is defined as 

d s (C)= min d S (V,V). (5) 

V,V'eC:VjtV 

A subspace code with minimum distance ds(C) > 2{p + t) is capable of correcting any p erasures and 
t errors with the minimum distance decoder. That is, if 2d$(V, U) < dg{C), the transmitted V can be 
recovered from the received U. 

One major construction of subspace codes [5J is through lifting the so called rank-metric codes 0. 
Let F™ xm be the set of all n x m matrices over F q . For 1,7 e F q xm , the rank distance between X 
and Y is defined as 

d R (X, Y) = rankpT - Y). (6) 

A rank-matric code M. is defined to be a nonempty subset of F™ xm . Each codeword of M. is a n x m 
matrix over F q . The minimum (rank) distance of M is defined as 

The most well-known rank-metric codes are Gabidulin codes @, which have the maximum possible 
minimum rank-distance, analogous to the Reed-Solomon codes in Hamming metric. 

Let I n be the n x n identity matrix. Denote by (X) the vector space spanned by rows of a matrix X 
over F q . The lifting of a rank-metric code M. gives the subspace code 



C = (v : V 



, X e M \. (8) 



In X 

It can be proved that ds(C) = dji(M.) f5\ . If M is a Gabidulin code, two efficient decoding algorithms 
have been developed for the resulting subspace code H, (5J. Both algorithms are guaranteed to decode 
up to the error-correction capability of the subspace code. 
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III. Layered Subspace Codes 



A. Code description 

Let Mi C F^ xm (l = 1,2, 



, L) be L rank-metric codes. We define the overall subspace code as 



C = < 



V : V 



i ni o ••• x 1 



In 



x 2 



(9) 



• • • I nL X L 

By lifting the rank-metric codes, we can obtain L component subspace codes 

C, = {*5:F, = ([(),, I ni l2 X,]),X,eM} a = l,2---,L), (10) 

where 0/ 1 is the n; x fX)|=i n «) all-zero matrix and 0; 2 is the n; x n «) all-zero matrix. For 

decoding purpose, we will assume the rank-metric codes to be Gabidulin codes. 

Obviously, for any Vi € C{ and Vj € Cj (i / j), ViDVj = {0}. Therefore, we have 

Property 1: C = {V : V = Vi V 2 © • • • V L , Vi 6 Cj} . 

The property leads to a superposition coding scheme, which is depicted in Fig. 1 . The overall subspace 
code consists of L superimposed layers (each corresponding to a component code), and hence we have 
the name layered subspace code. 



Encoder 1 



Encoder 2 



Encoder L 



V, 



V, 




V 



Fig. 1. The overall subspace code as a superposition coding scheme. 



Based on the definition (9) and (10), we further have the following property. 
Property 2: For any two codewords of C, V = V\®- ■ -®Vl {Vi € Q) and V = V[@- 
V = V' if and only if Vi = V{ for all I. 
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For the minimum distances of C and Ci, the following property holds. For the proof, see ll6l . 
Property 3: d s {C) = min {ds(Ci), d s (C 2 ), ■ ■ ■ ,d s (C L )}. 

B. Decoding algorithm I 

Suppose that a codeword V € C was transmitted and the vector space U is now received. Corresponding 
to each V\ C V, we define a subspace U\ C [/ as follows. It consists of all vectors of U such that the 
elements at the coordinates {l, • • • , Yli=\ n *} u { S!=i ra i + 1> ' ' ' > S^=i n i} are zero - Now, we can 
describe decoding algorithm I as follows: 

1) Extract Ui from U; 

2) Use the decoder for Ci to recover Vi from Ui 0, Q. 

Note that given a set of vectors that span U, a basis of U\ can be extracted with the aid of Gauss-Jordan 
elimination. Based on Property 2, once V/ G for all I can be recovered, the transmitted codeword V G C 
can be determined as V = V\ ® ■ ■ ■ ® Vl- 



Extractor 1 




Decoder 1 


► 








Extractor 2 


Decoder 2 


► 



U r 



Extractor L 



Decoder L 



-v, 



Fig. 2. The schematic diagram of decoding algorithm I. 



An illustration of the decoding algorithm is given in Fig. 2. It is seen that a parallel implementation 
is allowed. Moreover, if a receiver is only interested in a particular V\, then he only needs to perform 
the corresponding layer in Fig. 2. 

We now focus on the error-correction ability of decoding algorithm I. We first need to introduce a 
general result. 

Lemma 1: Let A and B be two subspaces of W. Let A be a subspace of A. Then, 

dim(A) - dimLA nB)> dim(A') - dim(A' n B). (11) 
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Proof: Since A n B is a subspace of A, there exists a (not uniqueQ) subspace C of A such that 
A = (A n -B) © C. Similarly, there exists a (not unique) subspace D of A' such that A' = (A' C\B)®D. 
Therefore, we only need to prove dim(C) ^ dim(D). 

Assume that there exists an x ^ £ D (~) (An B). Then x G D and x e B. Since A' = (A' C\B)@D, 
we have x & A . So x & A C\ B, which together with x G D contradicts the fact (A n 5) n D = {0}. 
Therefore, Dn(AnB) = {0}. Since A is a subspace of A (by hypothesis) and A' = (Ar\B)®D, D is a 
subspace of A. So (An5)©D is a subspace of A = (Af\B)@C. Therefore, dim(C) ^ dim(L>), and the 
statement of the lemma follows. □ 
For decoding algorithm I, we have the following result. 
Theorem 2: d s (U, V) ^ d s {U h Vj) for all I. 
Proof: It is important to note that 

v l nu = v l nu l (12) 

and 

u l nv = u l nv l . (13) 

Then, based on Lemma 1, we have 

dim(V) - dim(V nU)^ dim(Vi) - dim(^ n U{) (14) 

and 

dim(?7) - dim([/ n V) > dim(^) - dim(E/i n V t ). (15) 

Summing up (14) and (15), we obtain 

dim(V) + dim(U) - 2dim(£7 n V) > dim(V z ) + dim(ETi) - 2dim(C/ / n V t ). (16) 

By the definition of subspace distance, ds(V, U) ^ ds(Vi,Ui). □ 
Combining Property 3 and Theorem 2, we have 

Corollary 3: If 2d s (V, U) < d s {C), then 2d s {V u U{) < d S (Ci) for all /. 

'For an arbitrary w € A n B, we have C = {c + to : c G C} / C such that A = (inB)©C'. Note that all such 
vector spaces are isomorphic to the quotient space A\(A PI B) 1111 . 
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The corollary indicates that decoding algorithm I is guaranteed to decode up to the error-correction 
capability of the overall code. 

Note that it may happen that 2ds(V, U) ^ ds{C) while 2ds{Vi, Ui) < ds(Ci). In this case, the decoding 
algorithm can decode beyond the error-correction capability of the overall code. We give an example to 
show this. 
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The vector is filtered out. 
(a) For decoding algorithm I 
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A recovered vector in V, and the filtered-out 
vector result in an error vector. 



(b) For decoding algorithm II 



Fig. 3. A transmitted codeword V £ C and the corresponding received U. 



Example 1: Let the overall code C be composed of two component codes C\ and C%, which have 
parameters n\ = 3, ri2 = 4, m = 4, ds{C\) = 6 and ds{C2) = 8. According to Property 3, we have 
d${C) = 6. Fig. 3 (a) gives a transmitted codeword V € C and the corresponding received vector space 



s 



U. Also shown in the figure are V x € C x , V2 € C2, and U x and U2 extracted from U. We note that in 
extracting U x and U2, the vector of U in the last row is filtered out. 

It is easily verified that d s (V, U) = 7+7-2 x 5 = 4, d s (V x , U x ) = 3 + 3-2 x 2 = 2, and d s (V 2 , U 2 ) = 
4 + 4-2x3 = 2. Consequently, we have 2d s (V, U) > d s (C), 2d s (V x , U x ) < d s (C x ), and 2d s {V 2 , U 2 ) < 
<is (£2)- Therefore, for the instantiated V and U, decoding algorithm I decodes beyond the error-correction 
capability of the overall code. □ 

C. Decoding algorithm II 

It is well-known that for conventional superposition coding, the (iterative) successive interference 
cancellation (SIC) decoding algorithm is usually adopted ifTOl . Viewing a layered subspace code as a 
superposition coding scheme, we develop a SIC-like decoding algorithm, which is shown in Fig. 4. With a 
slight abuse of notation, the '©' here denotes the sum of two vector spaces that do not necessarily intersect 
trivially. By taking the dashed arrows into account, we obtain an iterative version of the algorithm. When 
decoder I does not decode into a codeword of Ci (this can be checked by the decoder), we set Vi to be 
the zero subspace. So if there only occur erasures in the operator channel, the iterative version in general 
outperforms its non-iterative counterpart. 



Extractor 1 



u; 



Decoder 1 



-9 U 



Extractor L-1 



Decoder L-1 



Extractor L 



u: 



u- 



Decoder L 



V, 



Fig. 4. The schematic diagram of decoding algorithm II (dashed arrows for the iterative version). 



It should be pointed out that in |8], the authors have used the idea of SIC to decode their constructed 
subspace codes. However, they did not mention any iterative decoding. 
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On the error-correction ability of decoding algorithm II, we have the following result. 

Theorem 4: If 2d s {V, U) < d s (C), then d s (V, U) ^ d s (V, U + V L ) ^ • • • ^ d s (V, U + V L + ■ ■ ■ + Vi). 
Proof: We prove the theorem by induction on I. Since 2ds(V, U) < ds{C), from Corollary 3, we 
have Vi = V t C V. Therefore, A\m{V + U) = dim(V+ + VJ)) anddim(FnC/) < dim(Fn (U + V t )). 
Based on the definition of subspace distance, we have ds(V, U) ^ ds(V, U + Vl). 

Assume that d s (V, U) > d s {V, U + V L ) > ■ ■ ■ > d s (V, U + V L + • ■ ■ + Vt) holds. Since 2d s (V, U) < 

ds{C), 2ds(V, U + Vl-\ \-Vi) < d s {C). From Fig. 4, we see that U + V L -\ h V/ is the input to the 

extractor l—l. Based on Corollary 3, we have Vi-\ = V\-\ C V. Consequently, ds(V, U+Vl-\ \-Vi) ^ 

ds(V, U + Vl + • • • + Vi-i). Thereby, the proof is complete. □ 

From the proving process, we see that decoding algorithm II is guaranteed to decode up to the error- 
correction capability of the overall code. 

Like decoding algorithm I, decoding algorithm II also occasionally decodes beyond the error-correction 
capability of the overall code. We note that the two algorithms may correct different errors in this case. 
We show this through the following example. 





Fig. 5. Decoding for another received U. 
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Example 2: For the parameters given in Example 1, we now use decoding algorithm II. Since U' 2 = U2, 
V2 can be recovered. The resulting U + Vj> and U 1 is shown in Fig. 3 (b). It is easily obtained that 
ds(Vi, U[) = 3 + 4 — 2x2 = 3. Since 2ds{V]_, U{) = ds(C\), V\ cannot be recovered with decoding 
algorithm II. 

Consider again the transmitted V in Example 1, but now suppose that the received U is given as in 
Fig. 5. It is easily verified that ds(V, U) = 7 + 6- 2x5 = 3, d s {Vi, U x ) = 3 + 2 - 2 x 1 = 3, and 
ds(V2, U2) = 4 + 4 — 2x3 = 2. Therefore, with decoding algorithm I, only V2 can be recovered. 

From the recovered V2, we calculate U+V2 and U[ as in Fig. 5 (b). Since ds(Vi, U^) = 3+3—2x2 = 2, 
we have 2ds{V\, U[) < ds(C). So V\ can be recovered with decoding algorithm II. 

In summary, for the transmitted V and the received U in Fig. 3, only decoding algorithm I can correctly 
decode, while for the same V and another U as given in Fig. 5, only decoding algorithm II can correctly 
decode. □ 

IV. Applications 

In this section, we show that layered subspace codes can be applied in various scenarios for network 
coding. 



Source 




Receiver 1 Receiver T 
(a) Single-source multicast 



Source 1 Source L 




Network 




Receiver 1 Receiver T 
(b) Multi-source multicast 



Source 1 Source T 




Receiver 1 Receiver T 
(c) Multiple unicasts 



Fig. 6. Possible application scenarios of layered subspace codes. 
• Single-source multicast 
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In this scenario, a single source communicates its information over a network to a specified set 
of T receivers, as shown in Fig. 6 (a). The source is encoded with the overall subspace code and 
the basis vectors defining a codeword V € C in (9) are transmitted. To deal with network dynamics, 
the number of component codes can be adapted. Therefore, this leads to an adaptive transmission 
scheme. On the other hand, by using component codes with different error-correction capabilities, 
the coding scheme can be used for unequal protection transmission [9]]. 

• Multi-source multicast Q 

As shown in Fig. 6 (b), L sources transmit independent information over a network to a specified 
set of T receivers. Source I is encoded with component code I and the basis vectors defining a 
codeword V\ € C\ in (10) are transmitted. Based on the received vectors, each receiver tries to 
recover V\ for all I. 

• Multiple unicasts 

A unicast means that a single source communicates its information over a network to a single 
receiver. In the multiple unicasts scenario, the number of receivers is equal to that of sources, as 
shown in Fig. 6 (c), and receiver I only requests the information from source I. The coding scheme 
is the same as in the multi-source multicast scenario. Since receiver I only wishes to recover Vi, 
decoding algorithm I is preferred in this scenario. 

V. Conclusion 

We treated the layered subspace codes in Q as a superposition coding scheme and proposed two 
efficient decoding algorithms. Error-correction abilities of both algorithms are analyzed. As an error 
control scheme, layered subspace codes can be expected to find various applications for network coding. 
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